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We investigate the influence of the spatial extent of the proton magnetization and charge densities 
on the 25 hyperflne splitting in muonic hydrogen. The use of a non-perturbative relativistic Dirac 
approach leads to corrections of 15% to values obtained from the perturbative treatment encapsu- 
lated by the Zemach radius, which surpass the next-leading order contribution in the perturbation 
series by an order of magnitude. 
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Precise values of the hyperfine splitting of the energy 
levels of electronic and muonic hydrogen have long been 
known to be relevant to precision tests of quantum elec- 
trodynamics [l|, Interest in the 25" state of muonic 
hydrogen has been strongly stimulated by the recent pre- 
cision measurement of the Lamb shift transition energy 
between the 2P^^^ and '2,8^^'^^ states by Pohl et al. 
This experiment led to the stimulating conclusion that 
the proton root mean square radius differs from the pre- 
viously accepted value in the literature by 4.9 a. De- 
termining the proton radius depended on extracting the 
Lamb shift of interest from the energy of the measured 

states. This 



transition between the 28^^^^ 



and 2P3^/=2 



extraction relied on a variety of mainly perturbative cal- 
culations performed by many authors, including those 
of MartynenkoJi[|L Borie @-[3, Pachucki @, Ofl, and 
many others 0, [Til- 13 • Our focus here is on the hyperfine 
splitting between the 25" states with total spin F = 1 and 
P = 0. A computed value (22.8148(78) meV 0) was used 
in . As a guide to the importance of precision, we note 
that a change in this value of 1% would correspond to a 
shift by 1(7 in the proton radius. Furthermore, this tran- 
sition is under current experimental investigation, with 
an announcement of a result expected soon. 

The F = and F — \ levels of atomic 8 states are split 
by the interaction between the magnetic dipole moments 
of the lepton and proton. Fermi computed this splitting 
treating the lepton and proton as point-like particles hav- 
ing magnetic moments but no spatial extent, and using 
the Coulomb potential to compute the non-relativistic 
hydrogen wave function. This first-order treatment gives 
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the Fermi energy of the nS states: 



E 



3n3 (m^ + Mp)3 ' 



(1) 



where = 2.792847351 is the proton magnetic moment, 
a is the fine structure constant, and Mp, to,, denote val- 
ues of the proton and muon masses. Zemach [l3| included 
the influence of the spatial extent of the proton charge 
and magnetization distributions using first-order pertur- 
bation theory. The relevant spin-spin interaction is given 

by 



Zemach 



12m ^Mp 



(2) 



where the scalar product of the lepton and proton Pauli 
matrices, cti • (T2 is +1 in the F = 1 state and -3 in 
the F = state, and pM is the magnetization density 
normalized to unity. A contribution of the muon anoma- 
lous magnetic moment is treated as a further (0.1%) cor- 
rection as per Ref. Q to which we compare. Zemach 
treated the difference between using Vzemach and the cor- 
responding interaction obtained by taking the magneti- 
zation density to be a delta function at the origin. His 
result can be expressed as 



2fia f dp 



^4 
P 



GE{~f)GM{-p^) 

Pp 



- 1 



EF{~2fia)Rp . 
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where Rp denotes the appropriately-named Zemach ra- 
dius of approximate value Rp = 1.040(16) fm Q and ji 
denotes the reduced muon mass. The appearance of 
results from its influence on the bound-state wave func- 
tion at the origin and Gm enters directly in the interac- 
tion. Effects of the spatial extent of the magnetization 



2 



density and higher order effects of Gm 7^ 1 arc not in- 
cluded. 

Our aim here is to treat the effects of Vzomach exactly 
within the framework of very precise numerical solutions 
of the two separate F = and F = 1 Dirae equations. 
The dominant binding interaction is taken to be the sum 
of the Coulomb potential and its lowest-order correction 
arising from vacuum polarization, each being modified by 
the non-zero spatial extent of the proton [15| . Since the 
size of the hyperfine splitting is much less than the size 
of the Lamb shift, a high degree of numerical accuracy is 
required. The method we use to do this is detailed in a 
previous publication |l5|. These high-precision calcula- 
tions have been shown |15| to reproduce exact, analytic 
eigenvalues for a point-Coulomb potential and to satisfy 
a virial theorem test to within 0.5 /icV which we set as 
an upper bound on our numerical errors here. 

To proceed we need to specify the magnetization den- 
sity pm and corresponding magnetic radius (r*^)^ = 
j (Pr r'^pM{r) (while including the effects of the charge 



radius (r^) 



C\2 — 



rp)c in the dominant binding poten- 



tials). The magnetization density is taken to be a nor- 
malized exponential 



[13 value = 0.8768 fm (re- 
cently updated, as per NIST) and that found in the anal- 



r^— the 2006 CODATA 



ysis of Pohl et al. 



0.84184 fm. This yields 



^EfP^A + B{{rl)M) + C{{rl)M?'\ 



(6) 



for which the values of A, B, and C are given in Table |T] 
for parameterizations at the two values of r^. 

The two parameterizations are compared in Fig. [1] in 
which we note that the dependence on r^^ is almost non- 
existent, while the dependence on is stronger. Also 
shown in Fig. [T] are the magnetic radii from Refs. [l8l.[T9| 
corresponding to the extractions from the Mainz (0.777 
±0.029 fm) and Jefferson Laboratory (0.850 ± 0.030 fm) 
collaborations, respectively. 

The large difference between extracted magnetic radii 
does not heavily influence the HFS. The difference be- 
tween using the Mainz value and the JLab value is only 
0.35 /ieV for the smaller value of the charge radius and 
0.23 p,eV for the larger value. Therefore, we can average 
the HFS obtained from the two values of the magnetic 
radius to calculate the value of the 25* hyperfine shift for 
each of the value of charge radius. The results are 



PM{r) 



V 

8^' 



77 = ^12/ {rl)M 



(4) 



The astute reader will question whether this exponen- 
tial form is a suitable choice for the magnetization den- 
sity. While the dependence of the Lamb shift energy 
on the shape of the charge distribution has been fully in- 
vestigated in Ref. [l^ — and discounted — this conclusion 
docs not automatically extend to the magnetization den- 
sity, which appears un- integrated in Eq. ([2]), as opposed 
to the integrated appearance in the finite- Coulomb and 
finite vacuum polarization potentials We shall re- 

turn to this issue, but for now an exponential form shall 
be assumed. 

We calculate the converged eigenvalues for the 25" hy- 
perfine splitting (being the difference between the eigen- 
values for the F — and F = 1 hyperfine states) by 
numerically integrating the Dirac equation 



AE. 



HFS,Z 

'2S 



22.6384 meV if = 0.84184 fm 
22.6347 meV if = 0.8768 fm 



(7) 



This is the value of the 25" hyperfine shift obtained 
using a complete treatment of the spatial extent of the 
proton's electromagnetic distributions. We can deter- 
mine the importance of our complete treatment as com- 
pared to using only the Zcmach radius by computing the 
2S hyperfine shift in the presence of only the finite-size 
Coulomb and finite-size vacuum polarization potentials, 
which we also fit to a polynomial form to give 

AEfP = 22.8521 - 0.1795 {rl)c + 0.0739 {rlfj^ meV, 

_ (8) 
as per Ref. [I^l- For the two chosen values of r^' this 
leads to an uncorrected (with respect to the finite size of 



dr [FMr)) " 1^ -A25 + V ^ ) 1^^25(0 j 

(5) 

for which K2s = —1, M again denotes the reduced muon 
mass, and A25 denotes the 25 eigenvalue shifted by the 
reduced mass. In doing so we obtain the upper and 
lower components {G2s{r) and F2s{r), respectively) of 
the muon wave-function in response to the combination, 
V of the finite-size Coulomb potential, finite-size vacuum 
polarization potential, and potential given by Eq. ^ 
(separately for each value of cti • (72), for a variety of 
values of and r*^. 

We express the computed hyperfine splittings (HFS) as 
a polynomial function of r^^ for two significant choices of 



TABLE I: Coefficients of Eq. ((H)) parameterizing the 25 
HFS with magnetic radius, relevant to the Zemach correc- 
tion for two values of the rms charge radius, calculated using 
the finite-Coulomb (Vfc), finite vacuum polarization (Vfvp), 
and magnetization potentials (Vzomach); and alternatively the 
finite-Coulomb potential plus magnetization potential. 



A [meV] B [meVfrn'^] C [meVfm-^] 

0.84184" 22.6085 0.1425 -0.1191 

0.8768° 22.6037 0.1450 -0.1203 

0.8768'' 22.5318 0.1388 -0.1142 

"Calculated using V = Vpc -|- Vpyp + Vzomach 
''Calculated using V = Vpc + Vzemach 
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FIG. 1: Comparison of data and parameterizations of Zemach corrections to the 2S hyperfine state of muonic hydrogen for 
two values of the proton rms charge-radius. Also shown are the calculated fits ('Fit 1' and 'Fit 2' denote fits for = 
0.84184 fm and = 0.8768 fm respectively) to the data as presented in the first two lines of Tableland the magnetic radii 
extracted from Bernauer et al. [11] [Mainz, {rliy^^ = 0.777 ± 0.029 fm] and Gilad et al. ,19] [Jefferson Lab HaU A GoUaboration, 
{r%j)^^^ = 0.850 ± 0.030 fm] used to calculate the Zemach corrections of Eq. ((7|. The point denoted by a star corresponds to 
the Zemach correction calculated via charge and magnetization densities extracted from Ref. ^20] at = 0.878 fm and 
r*^ = 0.860 fm as per Eq. (fTf)). 



the proton magnetization distribution) HFS of 

f 22.7690 meV if ?f = 0.84184 fm 
AS^H/s^J " . (9) 

[ 22.7639 meV if ?f = 0.8768 fm 

The isolated Zemach correction calculated here is given 
by the difference between the corrected calculation and 
the uncorrected one, 



Zemach 



= AE^^^^^ - AEfP 



We may compare this to the uncorrected shift calcu- 
lated using only the finite- Coulomb potential in Ref. [isj . 
interpolated to r?^ = 0.8768 fm, and we obtain 



AE^P 



22.6910 mcV, 



(12) 



indicating a correction using the finite- Coulomb potential 
alone (neglecting the vacuum polarization) of 



A£;zcmach = -0.1289 meV. 



(13) 



-0.1306 meV if 



0.84184 fm 



-0.1292 meV if rS, 



0.8768 fm 



(10) 



In order to identify the significance of vacuum po- 
larization on this result, we furthermore recalculate the 
Zemach correction using only the finite-size Coulomb po- 
tential plus potential given by Eq. ([2]). In this case, cal- 
culating the potentials using = 0.8768 fm for various 
values of r^^ we obtain a polynomial of the form given by 
Eq. ^ with values of A, B, and C listed in the last line 
of Table HI Using this polynomial form we calculate the 
HFS contribution for the two values of r^^ and average 
(though they are separated by only twice our numerical 
error) to obtain 



AE. 



HFS.Z 
2S 



22.5621 mcV. 



The deviation of this value from that in Eq. (fTO| for the 
same value of is a mere 0.3 fJ,eV, indicating that vac- 
uum polarization has negligible influence on the Zemach 
correction within our numerical error limits. 

Returning to the issue of the significance of the shape 
of the various distributions used in these calculations, 
we have recalculated the Zemach correction using charge 
and magnetization densities extracted from fits to form- 
factor data [l^l explicitly. The so-called charge and mag- 
netization densities pc,m('') may be extracted from the 
Sachs electric and magnetic form factors Ge,m{Q'^) via 
a Fourier transform 



PC, Mir) 



(27r)3 



''GE,Miq ) , 



(14) 



where for an atomic system we have = q , and we 
(11) take the Sachs form factors as the fits given in Ref. [l^j 
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of the form 



1 



1 + + QiT^ 



(15) 



1 +peT +pioT^ +Pi4r^ 



1 + P2T + PiT^ 



-P12T^ +Pl6T° 



(16) 

(noting that the proton magnetic moment /ip appears in 
Eq. ^ here rather than in Gm) with coefficients qi, and 
Pi given in Table HIl and for which t = Q'^/AM^. Us- 
ing the extracted charge density in the finite- Coulomb 
potential (which as per Ref. [l^l is constrained such that 
Tp = 0.878 fm), and the extracted magnetization density 
(constrained such that r^^ = 0.860 fm) in the Zcmach po- 
tential (the finite vacuum polarization potential, which as 
noted before is of minimal consequence, remains calcu- 
lated using an exponential charge distribution character- 
istic of ~ 0.8768 fm), we once again recalculate the 
Zemach correction to the 25* HFS, and we find 



AE^g^'^ = 22.6335 meV. 



(17) 



M 



As we do not have knowledge of either the or 
polynomial dependence of this energy, we may only com- 
pare to an interpolation of our previous results at the 
appropriate values. Comparing Eq. (|17p to our polyno- 
mial fit corresponding to = 0.8768 fm (the difference 
in rp being minimal) as per the second line of Table U 
and interpolating to r^^ = 0.860 fm, we find a difference 
of a mere 1.0 p-eY, indicating that the exponential distri- 
bution is indeed a good approximation for the magneti- 
zation density. This Zemach corrected HFS is included 
in Figure [1] where it is clear that it produces good agree- 
ment with our usage of a simplified model of the charge 
and magnetization densities. 

We compare our values of the Zemach correction 
(Eq. (fTO)) and Eq. (jlSp) to that which was used in the 



TABLE II: Coefficients of polynomial fit to Sachs electric and 
magnetic form factor data for the proton taken from 20j as 
used in Eqs. (fT5l[T6)) . 



Pi 



2 
4 
6 

8 
10 
12 
14 
16 



9.70703681 
0.00037357 

-1.43573 
0.00000006 
1.19052066 
9.9527277 
0.25455841 

12.7977739 



of 



TABLE III: Comparison of Zemach radii (Eq. (O) for 
and values of r*^ taken from 
[lil {r^' = 0.850 fm), and ^ 



several choices 
Refs. [il (r^* = 0.777'^ fm) 
{r^' = 0.860 fm). 



' p 


Rp (fm) [18] 


Rp (fm) m 


Rp (fm) ]20} 


0.84184 


1.022 


1.068 




0.8768 


1.046 


1.091 




0.878 






1.081 



analysis of Pohl et al. viz. that of Martynenko Q 
(denoted in that reference as 'Proton structure correc- 
tions of order a^') 



AE, 



Martynenko 
Zemach 



= -0.1518 meV. 



(18) 



14.5187212 

40.88333 
2.90966 

99.999998 

-1.11542229 
0.00004579 
0.003866171 

10.3580447 



This includes the influence of the Zemach term as well as 
the small (~ 5 %, depending on the value of the Zemach 
radius) effects of a putative treatment of the two-photon 
exchange interaction (with intermediate nucleon states). 

The comparison between Eq. (jlOp and Eq. ([T8|) sug- 
gests that the value used in the analysis of Pohl et al. re- 
quires a significant modification of some 10-15%. A more 
precise consideration of the importance of using a com- 
plete non-perturbative treatment is to look at the the 
next-leading contribution of Ref. [1] ('Proton structure 
correction of order a^' = -0.0017 meV). This is an order 
of magnitude smaller than the difference found here, indi- 
cating that the perturbative approach is not satisfactory 
in this scenario. 

We may also study the importance of the Zemach ra- 
dius by examining the values of the Zemach radii for the 
values of r^'*^ of present interest, as displayed in Ta- 
ble mil The HFS is almost independent of the value of 
Tp^, but the Zemach radius displays a very significant 
dependence on this parameter. 

The difference between our value quoted in Eq. (fTUl) 
and the value of Martynenko quoted in Eq. (fT8|) docs not 
resolve the proton radius puzzle. The effect presented 
here alters the 25* hyperfine splitting by 0.02 meV, con- 
tributing only 0.005 meV (while 0.3 meV is needed) to 
the transition measured by Ref. 0]. However, our results 
add weight to the argument that the perturbative cor- 
rections require further attention and re-evaluation using 
the latest numerical tools, along with a close examination 
of the relevant physics to ensure that no contribution is 
overlooked, such as that reported in Ref. [2]| . 

Recomputing the 25* HFS with this new Zemach in- 
clusion, we obtain a splitting as given in Table IIVI which 
summarizes our updated non-perturbative calculation of 
this splitting. 

We note a relevant error in Ref. 
tributions to the 25 HFS labelled 'Proton structure cor- 
rections of order and a^' are listed in Table HI of 
Ref. [1^ as included by the Dirac calculation, whereas 



151 — the Zemach con- 
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these terms should appear under 'Remaining corrections'. 
This contribution now appears corrected in Table IIVI 
which should serve as a replacement to the former ta- 
ble. The re-calculation of the Zemach contribution alters 
the prediction of the proton rms radius via the analysis 
of Ref . [m to arrive at 

0.84182(67) fm, (19) 

which does not differ from the former conclusion 
{rp = 0.83811(67) fm) in a statistically significant fash- 
ion with regard to the main discrepancy, but does agree 
with that of Pohl et al. to a high level of precision. 

We conclude by stating that our calculation of the in- 
fluence of the spatial extent of the proton's electromag- 
netic distributions using a non-perturbative, relativistic 
framework results in a 10-15% deviation of the Zemach 
correction from the value used by Pohl et al.. Further- 
more, the very weak dependence of the HFS on the value 
of the magnetic radius, shows that the perturbative treat- 
ment, dominated by the influence of the Zemach radius, 
does not capture the effects of the proton's spatial distri- 



butions on the HFS. This result invites and encourages 
further investigation into the many corrections which en- 
ter the analysis of the muonic hydrogen spectrum at the 
level of precision required to extract the value of the pro- 
ton radius at the required level of precision. 
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TABLE IV: Contributions to the 2Si/2 hyperfine splitting with comparison to values found in Martynenko 01 corrected and 
updated from Ref. 15] to include the re-calculated Zemach contribution. Values are all in meV. Errors in the Dirac calculations 
are taken to be ±500 neV. The listed corrections are already included in our Dirac calculations and are listed by their descriptions 
in Ref. AH further corrections to both the perturbative calculation and our calculation are contained in 'Remaining 

Corrections' which encompasses the muon AMM, amongst other corrections listed in Ref. Q. We note that the in Ref. Q 
the Zemach correction is listed as 'Proton structure corrections of 0{a^)' and may not include considerations of finite-size 
in the wave-function. The polynomial dependence on (r^)" of this splitting is not discussed in the literature. The strict Vp 
dependence of the Zemach correction has not been calculated, and as such we have selected a value (that corresponding to 
r-p — 0.8768 fm) suitable to our previous analysis 15]. 



Contribution 


Martynenko 




Present \Vork 




Dirac (V = Vc) 






22.8229 




Dirac (V = Vc -f Vvp) 






22.8976 




Dirac (V = V"fc) 




22.7774 


- 0.1746 (rp) + 0.0709 




Dirac {V = Vfc + Vvp) 




22.8510 


- 0.1701 (r^) + 0.0667 




Dirac {V = Vfc + Vfvp) 




22.8521 


- 0.1795 (rl) + 0.0739 




Fermi Energy Ef 


22.8054 








Relativistic correction ^{Za)'^ Ep 


0.0026 








VP corrections of orders a^, a® 










in the second order of perturbation series 


0.0746 








Proton structure corrections of order 


-0.1518 








Proton structure corrections of order 


-0.0017 








Zemach correction 






-0.1292 




Subtotal: 


22.7291 


22.7229 


- 0.1795 (r^) + 0.0739 




Remaining Corrections 


0.0857 




0.0857 




Total: 


22.8148 ± 0.0078 


22.8086 


- 0.1795 (rl) + 0.0739 





